CP Non-conservation in the leptonic sector 



Petre Di^ 

In this paper we use an exact method to impose unitarity on moduli of neutrino PMNS matrix 
recently determined, and show how one could obtain information on CP non-conservation from a 
limited experimental information. One suggests a novel type of global fit by expressing all theoretical 
quantities in terms of convention independent parameters: the Jarlskog invariant J and the moduli 
\U a i\, able to resolve the positivity problem of |{/ e 3|- In this way the fit will directly provide a value 
for J, and if it is different from zero it will prove the existence of CP violation in the available 
experimental data. If the best fit result, \U e 3\ 2 < 0, from M. Maltoni et al, New J.Phys. 6 (2004) 
122 is confirmed, it will imply a new physics in the leptonic sector. 



PACS numbers: 12.15 Hh, 12.15 Ff, 13.30 Ce, 23.40 Bw 
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I. INTRODUCTION 

Neutrino experiments have shown that neutrinos have 
mass and oscillate, the last property suggesting a mixing 
of leptons similar to that of quarks, Q-B; for more ref- 
erences on experimental data see, e.g., ||| and The 
presence of the lepton mixing opens the possibility to see 
a CP violation in the leptonic sector, phenomenon which 
is considered an indispensable ingredient that could gen- 
erate the excess of the baryon number of Universe 0, 
since the baryon number can be generated by leptogene- 
sis 8J. Excluding the LSND anomaly, all neutrino data 
are explained by three flavor neutrino oscillations, and 
the determination of oscillation parameters is one of the 
main goals of phcnomcnological analyzes. Recently two 
detailed combined analyzes of all neutrino data have been 
published, see Refs. and [||, that materialized in 
upper and lower bounds on sin#jj, these bounds being 
converted in Ref. Q into intervals for the moduli of the 
PMNS unitary matrix [l(J • In both above papers the CP 
phase 5 was not considered a free parameter, e.g., only 
the particular cases cos<5 = ±1 have been considered in 
la . We shall write the numerical moduli matrix as 



V 



0.835 ±0.045 
0.355 ±0.165 
0.365 ±0.165 



0.54 ±0.07 
0.575 ±0.155 
0.590 ±0.15 



< 0.2 

0.7±0.12 ) (1) 
0.685 ±0.125 



where the central values are given by half the sums of 
lower and upper bounds entering and the errors cor- 
respond to a three sigma level, excepting V e s for which 
only an upper bound is given. Complementary infor- 
mation concerns the absolute value of the two neutrino 
mass-squared differences |5j , but we have no information 
on: 

• the magnitude of C/ e 3-element of the leptonic mixing 
matrix, or the value of the mixing angle #13; 

• the existence of leptonic CP violation, or a value for 



the Dirac phase 5 7^ 0, n entering the PMNS matrix [To| : 

• the order of the mass spectrum, or the sign of Amf 3 ; 

• the magnitude of neutrino masses, etc.; see, e.g., [Til ]- 


Although it is believed that only the next generation 
experiments will give an answer to the first above items, 
|l4|. we consider that by using unitarity we can obtain 
some information from experimental data as given in Q 
for both the existence of CP violation in leptonic sector, 
and the parameter (9 13 , all the more as measurement of 
#13 is considered a first mandatory ingredient for inves- 
tigation of CP leptonic violation in — > v e transitions 
and for the mass hierarchy determination. Concerning 
CP violating phase, 6, the general opinion is that it will 
require a major experimental effort because of the intrin- 
sic difficulty to measure it, see e.g. Ref. |l5j . 

The current phenomenological analyzes of data are 
given in terms of mixing angles 0y, see, e.g., Q and [||, 
and phase 5 which all depend on the chosen parametriza- 
tion of the unitary matrix one starts with. In the pa- 
per we propose another method for analyzing the data 
by expressing all measurable quantities directly related 
to oscillation amplitudes in terms of moduli, |J7 a j|, and 
Jarlskog invariant J, fl6j ]. From a theoretical point of 
view this approach has the advantage that the recon- 
struction of a unitary matrix from knowledge of moduli 
leads essentially to a unique solution, see, e.g., Ref. [l7| . 
From an experimental point of view the method allows 
to get from nowadays limited information on neutrino 
physics all moduli of the unitary matrix such that one 
can predict the approximate form for the oscillation am- 
plitudes, for example v T — > f M , which could be of great 
help for experimenters. It is easily seen that the experi- 
menters have measured quantities directly related to os- 
cillation amplitudes, which contain information on the 
moduli \U e i\, I £42 1, I j , and | C7%- 3 1 . The above mod- 
uli are correlated by unitarity to | /7 e3 \ 



\U, 



e3 



= 1 — 117: 



el 



-117, 



<:'2 



= 1 - 
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\U^-\U T tf (2) 

By using J5J, from the first row and from the third col- 
umn in Q one gets the bounds 

- 0.1465 < |[/ e3 | 2 < 0.155, -0.3285 < \U e3 \ 2 < 0.35(3) 
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The moduli of lower and upper bounds in the above 
relations are of the same order of magnitude, which 
means that the possibility V* 3 < is not excluded. In 
fact the best fit of the KamLAND data in Ref. Q gave 
V e 3 = —0.13; see also Figure A. 2 in Ref. [j| and, respec- 
tively Figure 14 in Ref. ,6] . Hence the phenomenologists 
have to take into account this possibility, which, if con- 
firmed, could be interpreted as a breaking of the estab- 
lished rules, signaling really a new physics. A slight mod- 
ification of the theory can take it properly into account. 
After non vanishing masses, such a surprise coming from 
neutrino physics is not a priori excluded. 

As we will show in paper our approach makes full use 
of unitarity property, which implies many correlations 
between all moduli of a unitary matrix. For a complete 
determination of an (exact) unitary matrix it is neces- 
sary the knowledge of four independent moduli, see, e.g., 
Ref. Taking into account that a 3 x 3 unitary ma- 
trix has nine moduli, the choice of four ones is completely 
at our disposal. Hence for analyzes we may use moduli 
which, apparently, do not appear in the measured os- 
cillation amplitudes. In fact there are 58 groups of in- 
dependent moduli, and this fact, coupled with unitarity 
constraints, gives us an opportunity to find all moduli, 
even nowadays we have only a limited experimental in- 
formation. 

In the present paper we explore the possibility to ex- 
tract information on CP violation in the leptonic sector 
starting with the nowadays knowledge on moduli matrix 
(fj"|). We find that within the moduli intervals given by 
||TJ there are values corridors that are compatible to the 
existence of at least one unitary matrix, in fact one can 
construct a continuum of solutions, result which could be 
interpreted as a signal of CP violation in the leptonic sec- 
tor. A clear signal of CP violation will be a value J ^ 
for the Jarlskog invariant, value which can be provided 
by a global fit suggested in the paper. If the new fit will 
confirm the result, | t/e3 1 2 < 0, by Maltoni et al, 0, this 
will be equivalent with the existence of new physics in 
the leptonic sector. 

The theoretical tool that will be used in the following 
is unitarity of the PMNS matrix, and the constraints 
on physical observables generated by this property. We 
make full use of phase invariance property to write the 
matrix under the form 

U = 

C12C13 C13S12 S13 \ 

— C'23Si2e tS — C12S23S13 Cl2C23e J<5 — S12S23S13 S23C13 (4) 
Sl2S23,e lS — C12C23S13 — Cl2S23e l<5 — S12C23S13 C23C13 / 

and the notation is Cy = cos dij and Sij = sin dij , where 
Oij are the mixing angles with ij = 12, 13, 23, and 6 is 
the Dirac phase which encodes CP violation. 

The matrix U is assumed and built as a unitary matrix, 
but nobody assures that the moduli \Uij\, the angles %, 
or the phase 5 extracted from experiments come from a 
unitary matrix. Usually it is assumed that irrespective 
of how the measured data are, they are compatible to 
the existence of a unitary matrix. A problem which was 



not considered in physical literature is the property of 
unitary matrices to be naturally embedded in a larger 
class of matrices, the set of double stochastic matrices. 

A 3 x 3 matrix M is said to be double stochastic if its 
entries satisfy the relations [2(j 

3 3 

rriij > 0, ^ m i] = !) X! mi i = = 1 ' 2 ' 3 ( 5 ) 

»=i j=i 

All such matrices form a convex set, i.e. if Mi and M2 
are double stochastic, matrices from the set 

M 3 = x Mi + (1 - x) M 2 , with < x < 1 (6) 

are double stochastic [2(j. The unitary matrices are a 
subset of this larger set if we define the entries rriij of the 
double stochastic matrix by relations 

ma = Pitf (7) 

One easily sees that relations (jSJ are satisfied because of 
the unitarity property 

UU* = U^U = h (8) 

where f denotes the hermitian conjugate matrix, and I3 
is the 3-dimensional unit matrix. The subset of dou- 
ble stochastic matrices coming from unitary matrices is 
known as unistochastic matrices. It is well known in the 
mathematical literature that there are double stochastic 
matrices that do not come from unitary matrices, see, 
e.g., [2(j - 0. This unpleasant feature rises a few prob- 
lems which have to be solved before trying to do a fit of 
the experimental data. The first one is to find the neces- 
sary and sufficient conditions the data have to satisfy in 
order they should come from a unitary matrix. In other 
words we have to find the separation criteria between the 
two sets. After that we have to provide an algorithm for 
the reconstruction of U from the given data. With these 
aims in view we have developed a formalism for express- 
ing the necessary and sufficient conditions for a 3 x 3 
double stochastic matrix to come from a unitary one. 
On the other hand the embedding relation J7J suggests 
that a "natural" parametrization of 3 x 3 unitary matrix 
could be that given by four independent moduli. This 
can be done with a proviso: cos<5 as function of moduli 
must take physical values, i.e. — 1 < cob 8 < 1. The use 
of moduli is also justified by the fact that mixing angles 
and CP phase are not well defined objects, they depend 
on the form of the unitary matrix. In this respect there 
are at least two popular forms: the PDG proposal [24| . 
and the classical proposal [25| by Kobayashi-Maskawa, 
and the physics has to be invariant, i.e. not dependent of 
one or another choice. The moduli are such invariants. 
The above statements can be formalized by the follow- 
ing (obvious) axiom: irrespective of the physical processes 
where the moduli appear, or the phenomenological models 
used in analyzing the data, the numerical moduli values 
should be the same. 
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We used the novel formalism to develop an algorithm 
for the reconstruction of a unitary matrix from a double 
stochastic one when it is compatible to a unistochastic 
one. In this sense we modified the standard x 2- test to 
allow the implementation of all unitarity constraints, for 
a better processing of experimental error affected data. 

To see that the found unitarity constraints are nec- 
essary to be used in any phenomenological analysis we 
constructed an exact double stochastic matrix, from the 
moduli , which is perfectly acceptable from an experi- 
mental point of view, i.e. all the moduli are very close to 
the central values in being within the experimental 
error bars, but which is not compatible with a unitary 
matrix. 

Last but not least we used the convexity property of 
double stochastic matrices to develop a method for do- 
ing statistics on unitary matrices through their squared 
moduli, problem which was an open one in the literature 
until now, see, e.g., Refs. (2||, (231 ■ 

The paper has the following structure. In section 2 we 
present our phenomenological model that makes a novel 
use of unitarity property of PMNS matrix, the main point 
being the focus on the condition that separates double 
stochastic matrices from unitary ones. This condition 
is, in the same time, the consistency condition between 
the theoretical model, represented by the PMNS matrix, 
and the experimental data, as those appearing in l(T]l. 
In section 3 we present the reconstruction algorithm of 
unitary matrices from experimental data, when they are 
compatible, and show how errors affected data can be 
consistently included in the formalism. We propose a 
two terms x 2_ test which has to contain a piece enforcing 
the fulfillment of unitarity constraints, and another piece 
which will properly take into account experimental data. 
In section 4 we propose a new type of global fit, done in 
terms of invariant parameters, J and moduli |£/y|, which 
can resolve the existence of CP non-conservation in the 
lcptonic sector by using the nowadays experimental data, 
or, alternatively, prove the existence of new physics in the 
lcptonic sector. In particular we show that around central 
values from data Q there is a continuous (approximate) 
unitary set of matrices all of them being consistent with 
5 w 7r/2, for a broad range of numerical V e3 values. By 
changing the moduli values around one looks for unitary 
matrices one gets other values for S and V e z. The paper 
ends by Conclusions. 



II. PHENOMENOLOGICAL MODEL 

By phenomenological model we will understand in the 
following a relationship between the entries of a double 
stochastic matrix, provided in general by experiments, 
and the entries of a unitary matrix. The experimental 
data will be the numbers affected by errors entering QJ. 
We remind that in order to find allowed ranges for the 
moduli \Uij\, as those given in l[T)l. the solar neutrino and 
atmospheric neutrino measurements, as well as the reac- 



tor neutrino experiments have been used, see [f| , |(| , and 
The analyzes lead to numbers for all moduli, but 
\U e s\, from the first row and the last column. Unfortu- 
nately only three moduli are independent, because the 
unitarity gives 



el 



UelV = \U, 



H3 



\u t3 



(9) 



Taking into account the embedding Q, which does not 
imply the unitarity fulfillment by experimental data, we 
will use different notation for the theoretical quantities, 
and for experimental data. Hence in the following we as- 
sume that the data are contained in a numerical matrix, 
V, whose entries are positive 



V 



fV el V e2 V e3 
V 

K v Tl v t2 v t3 



(10) 



and the nowadays experimental data on the moduli 
V a i, a = e, /i, r, i = 1,2,3, are summarized in matrix 
Q). From a phenomenological point of view the main 
problem to be solved is to see if from a numerical matrix 
as Q one can reconstruct a unitary matrix as Q. For 
that we define a phenomenological model as follows 



V 2 = |£/f 



(11) 



relation that has to be understood as working entry-wise 
leading to the following equations 



V-2 2 2 T/"^ ^ ^ T/^ 2 

el — c 12 c 13) V e2 — s 12 c 13j V e3 ~ S 13 



fi3 



2 2 T/-2 _ 2 2 
s 23 c 13j v t3 — c 13 c 23i 



V, 



/'I 



S12C23 + Sl3 s 23 c 12 + 2S12S13S23C12C23 COS 8, 



V, 



M 2 = C12C23 + S12S13S23 - 2S12S13S23C12C23 cos<5, (12) 



VtI = Sl3 c 12C23 + S12S23 — 2S12S13S23C12C23 COS 5, 

2 2 2 2 2 2 

Vt2 = S12S13C23 + C 12 S23 + 2S12S13S23C12C23 COS (5 

together with the double stochasticity relations 

|2 -1 



i=l,2,3 

E \u a3 



£ Vai - 1 = 0, a = e, /J,, v 

i=l,2,3 

£ V&-1 = 0, J = 1,2,3 (13) 



The above phenomenological model is similar to that 
proposed by Wolfenstein, [2q|. i.e. it is a direct relation- 
ship between the measured values Vij and the theoreti- 
cal parameters Sjj and 5 entering U. The difference be- 
tween the two approaches is the explicit use of the double 
stochasticity properties (|13|l , and making no approxima- 
tions on right hand side equations (|12fl . Also important is 
the fact that on left hand side of equations (|12|l we have 
sets of nine numbers, Vij, which all could be obtained 
from experiments, and on right hand side there are four 
independent parameters. Hence checking the consistency 
of the system of equations Ijl2(l is a natural problem and 
it has to be resolved. Since the matrix U is parameterized 
by four independent parameters, e.g. Sy and S, there is 
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not at all obvious that the nine equations i|12|) have a 
physical solution for all experimental values allowed in 
i ffl l. On the other hand we have to be aware that the 
parameters sy and S are convention dependent and have 
no intrinsic theoretical significance, and this fact rise the 
problem of using only invariant functions generated by 
the entries of the matrix Such functions are the 

Jarlskog invariant J and the moduli V a i, see e.g. p9j . 
and this is the main reason for starting with moduli as 
independent parameters in our approach. 

We said before that for n > 3 there are double stochas- 
tic matrices which are not unistochastic, see, e.g. [2p|- 
|23|. hence a novel problem arise: what are the necessary 
and sufficient conditions for their separation. In this re- 
spect there are only two mathematical references, j2jj 
and |22j, where theorems are given on the necessary and 
sufficient conditions for a double stochastic matrix to be 
also unistochastic. See also Refs. 01 an d for the 
physicists point of view. The theorems are only existence 
theorems, i.e. they do not provide a constructive method 
for recovering a unitary matrix from a double stochastic 
one when they are compatible. However they are impor- 
tant since show the equivalence between unitarity and 
the positivity of (any) area generated by the so called 
unitarity triangles, which is one form of the consistency 
condition between double stochastic matrices and uni- 
tary ones. Here, by using the phenomenological model, 
Eas.^li p -(|13 p . we provide an alternative method for the 
separation criteria between the two sets, the consequence 
being a new form of the consistency conditions between 
the experimental data and the theoretical model, which 
leads eventually to a constructive algorithm for recovery 
of unitary matrices from error affected data. 

To see how the phenomenological model defined by the 
relations (fll|) - (|13|) works, let simplify the things and as- 
sume for a moment that the relations (|13l) are exactly 
satisfied by the experimental data, implying the data are 
not affected by errors. Then it is an easy matter to find 
from the first five relations <|12|l three independent ones 
which give a unique solution for the sy, ij = 12, 13,23. 
In other words, if the experimental numbers satisfy the 
relations 

v e \ + K 2 2 + v e % = i 

+ Kz + V T \ = 1 (14) 
we get always a physical solution for which is unique 



under the condition < s^- < 1, and depends on the 
three chosen independent moduli. Substituting this so- 
lution in the last equations one gets four equations for 
cos (5, that lead to a unique solution for it. But nobody 
guarantees that the solution will satisfy the physical con- 
straint 

-1<cos<5<1 (15) 

Thus the necessary and sufficient conditions that guaran- 
tee the consistency between the data and the theoretical 
model are given by: 

0<% <1, ij = 12, 13, 23, -1<cos<S<1 (16) 
The constraints i|16|) together with relations l|12|l - l|13|) 
prove the unitarity property of data, and provide an al- 
gorithm for recovery of unitary matrices from their ex- 
perimental moduli. 

The embedding relations Q of unistochastic matrices 
into the double stochastic set show that the number of 
independent parameters is the same for the two sets, i.e. 
in our case it is equal to four. Now we construct a double 
stochastic matrix by using the constraints l|13|) and the 
entries from the data set, (|ll)fl . In order to simplify the 
formulas form we make the notation 

V e i = a, V e2 = 6, V^i = c, V m2 = d, 

V^ = e, V Tl =g, V T2 = h, V T3 = f (17) 

As we already remarked the measured parameters a, 6, e, 
and / are not independent since unitarity implies the 
relation a 2 + b 2 = e 2 + f 2 . Thus in principle we have 
to chose only three of them. Our choice for the first two 
matrices is, a and 6, since they have the smallest errors, 
and c and e, respectively, d and /. The third matrix will 
be generated by c, d, g, h, i.e. by moduli that are not yet 
measured. If the independent parameters are a, 6, c, e, 
the matrix S 2 = (Sfj), where 

S = (Sij) = 

/ a b Vl - a 2 - b 2 \ 

c VI - c 2 - e 2 e (18) 

\ \/l - a 2 - c 2 \Jc 2 + e 2 - b 2 V 'a 2 + b 2 - e 2 J 

is double stochastic. Using the relations V 2 ^ = S 2 ^ in 
Eas. lfT^j) we find values for and cos£, e.g. 



cos 



S (l) 



= \J 1 — a 2 — b 2 , and 



Va 2 + b 2 1J " Va 2 + b 2 

-b 2 (a 2 + b 2 ) + c 2 (a 2 + b 2 ) 2 + e 2 (b 2 (a 2 + 1) + a 2 (a 2 - 1)) 



(19) 
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In the second case when the independent parameters are a, b, d, f one gets 



2 



,g> = -J = , a W = ^l-a 2 -^, and .sg) = V" 2 +^/ : 
12 Va 2 + b 2 V« 2 + & 2 

rf2l + & 2 )(& 2 (1 - a 2 - & 2 ) - d 2 (a 2 + b 2 )) + f 2 (a 2 + b 2 (a 2 + b 2 - 1)) 

cos<r ' = — — , 

2a6/Vl - a 2 - 6 2 ^/a 2 + b 2 - f 2 

and, similarly, in the third case when used parameters are c, d, g, h 



(3) Vl - d 2 - /l 2 (3) / 2 , , — 2 , , 2 r (3) Vl - c 2 - d 2 

12 V2 - c 2 - rf 2 - /i 2 13 V ^ 23 V2 - c 2 - d 2 - /i 2 



(20) 



(21) 



cos^ 3 ) = (c 2 (l - c 2 ) - d 2 (l - d 2 ) - g 2 {\ - g 2 ) + h 2 (l - h 2 ) + (d 2 g 2 - c 2 h 2 )(c 2 + d 2 + g 2 + h 2 - 2)) / (22) 
(2 Vl - c 2 - d 2 y/l - c? ~ g 2 yj\-d 2 - h 2 y/l - g 2 - h 2 V ' c 2 + d 2 + g 2 + h 2 - lj 



In the first two cases the mixing angles Sij depend only By using the form l|19J) for cos 8 the above condition is 
on three parameters, such that the only constraint on all equivalent to 
four parameters is given by — 1 < cos<5 < 1 which we 
write it in the form 

1 - cos 2 8 = sin 2 8 > 



• 2 c 

sm — 

-(a 2 + b 2 ) 2 {b 2 (b 2 -2c 2 (a 2 + b 2 )) + c 4 (a 2 + b 2 ) 2 + e 4 (l - a 2 ) 2 + 2e 2 b 2 (a 2 - 1) + 2c 2 e 2 (a 2 (a 2 - 1) + b 2 (l + a 2 )) 

4a 2 b 2 e 2 (l -a 2 - b 2 )(a 2 + b 2 - e 2 ) 

I 

On the other hand Branco and Lavoura, 01 > an d, re- 
spectively, Jarlskog and Stora, [19j . found a geometric 
invariant, the area of any unitarity triangle, as a measure 
of CP-violation. This area has to be a positive number 
if the CP symmetry is violated. If k, i=l,2,3, denote the 
lengths of the triangle sides, from l|18|l one gets 



> (23) 



e', 



h = a b, I2 = c yl 
l 3 = Vl-a 2 -c 2 Vc 2 + e 2 -b 2 (24) 



I 

One makes use of Heron's formula A = y/p{p — h)(p — h)(p — h), where p = (h + h + fa)/2 is the semiperimeter, 
for getting area A, and the reality condition A 2 > has the form 

16yl 2 = -(b 2 {b 2 -2c 2 (a 2 + b 2 )) + c(a + b 2 ) 2 + e 4 (l - a 2 ) 2 + 2e 2 b 2 (a 2 - 1) + 2c 2 e 2 {a{a - 1) + 6 2 (1 + a 2 ))) > 0(25) 

I 

Looking at equations (|23|l and 125(1 we see that they are condition (|16fl on cos 8 is a novel separation criterion be- 
equivalent, the supplementary factors appearing in i|23|) tween the unistochastic matrices and the double stochas- 
being positive for a double stochastic matrix. Hence the tic ones. The separation border is given equivalently by 
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A = 0, or cos 2 6 = 1. 

E quiv alent conditions to were given in 0, 
and , where positivity of sin 2 5 >0 was mainly used to 
constrain the four implied moduli. One aim of the paper 
is to put together all that information for devising a fit- 
ting model which will use all the possible constraints com- 
ing from unitarity. The idea behind it is to use the non 
uniqueness of the set of four independent moduli used to 
parameterize a unitary matrix which leads to many re- 
lationships between them and which over constrain the 
error affected experimental data. A consequence will be 
that unitarity implies a very precise tuning between all 
moduli, and we could get information on the physical 
relevant parameters even from a limited experimental in- 
formation. 

From (|19|) . i(20|) and ((21(1 we see that the expressions 
defining and 6 are quite different. The interesting 
thing is that all the possible forms for and cos 6 take 
the same numerical values when the numerical matrix is 
double stochastic. Let us consider the following matrix 



(26) 





( 


167 
200 


11 

20 


vTI 
200 


s = 




8 

25 


\/237 
25 


IS 
25 




\ 40 


/ 1603 

v — 


V3183 
100 


V19253 
200 



0.835 0.55 0.0166 
0.32 0.6158 0.72 
0.4476 0.5642 0.6938 

whose entries are within the error bars from ft}, and 
the matrix S 2 = (£?■) is a double stochastic matrix. 
The numbers entering it have been chosen such that the 
measured parameters, V e \, V e 2, and V T 3, should be 
around the central values from For this matrix we 
get 



„(0 



S 23 



110 



V39989 
144 



V39989 



0.55, 



0.72, 



'13 



200 
1,2,3 



0.0166, 



(27) 



5419565144 



2066625V211783 



-5.698, i = 1,2,3(28) 



The numerical results show that all the three angles Sy 
are physical, but S is not, i.e. the data matrix S, Eq. (l2"6l . 
does not come from a unitary matrix. Computing the 
triangle area we find 



A = 0.0107i 



(29) 



i.e. an imaginary number which send the same signal, 
namely the above data 1)26(1 are not compatible to the 
existence of a unitary matrix. 

Thus the new unitarity constraint — 1 < cos 6 < 1 puts 
non trivial constraints on the experimental data. From 
an experimental point of view the above matrix is per- 
fectly acceptable, from a theoretical point of view it is 



not; and the non-unitarity of the data can be checked by 
the above method. 

In case when the numerical matrix Q is double 
stochastic, as the above matrix S, s,y and cos<5 do not 
depend on the chosen four independent moduli used to 
define it. If the modulus of cos 5 is outside the inter- 
val (—1,1) the data are not compatible to the theoreti- 
cal model and the story ends here. If cos<5 takes values 
within the interval (—1,1), the data are compatible with 
the theoretical model, and the values for 6 and Sij nec- 
essary for the recovery of the unitary matrix are easily 
obtained, e.g. from the relations ((19(1 . and by introducing 
them in Eq. W e explicitly recover the unitary matrix. 

From almost any numerical experimental data as 
we can form a doubly stochastic one by choosing at our 
will four independent parameters, the other five being 
completely determined by equations similar to 1(18(1 . In 
this way we obtained the equations ((19(1 . 120|) and l(21|) 
for cos 6. By using numerical information from matrix 
(JTJ for a, b, c, d, e, g, h and /, we get 

cos^ 1 ' = -1.1 i, cos8 ( c 1} = -0.553, cos5 w = -0.985 
cos<5<_ 2) = 0.94 i, cos 5 ( c 2) = 0.118, cos 8^ = 0.737 (30) 
cos5^ 3) = 0.05 i, cos 5i 3) = 0.018 i, cos 8 (3) = 0.0027 i 

where + , c, — denote values of cos 5 function obtained 
from central values+cr, central values, and, respectively, 
central values— a. The corresponding values for cos<5 on 
the three rows arc different because the numerical val- 
ues for parameters a, b, c, e, respectively, d,f,g,h do 
not come from the same doubly stochastic matrix, e.g. 
e V 1 — c 2 — d 2 , or numerically 0.7 ^ 0.737. The above 
numerical computations send a clear signal: unitarity of 
the PMNS matrix is a property shared by all its entries 
implying strong correlations between the numerical val- 
ues taken by Vij , and this unitarity feature has to be well 
understood by people doing phenomenology. The border 
of the physical region, A = 0, or cos S = ±1, can be 
used to find maximal intervals for the unmeasured mod- 
uli. If we consider V e i, V e 2, and or V T 3 fixed to the 
central values one gets V^i € (0.31, 0.45) which represent 
40% from the 3 a interval given by (JIJ, and respectively, 
Vfj.2 S (0.53, 0.64), i.e. about 35% from the 3 a inter- 
val. If one computes now the same quantities by using 
the non unitary double stochastic matrix 1(26(1 one finds 
Vfii G (0.37, 0.39), and respectively, € (0.57, 0.59), 
i.e. intervals up to seven times shorter than in the pre- 
vious case. In conclusion the constraints are stronger for 
matrices which are far from unitary ones. The previous 
computations show also how one can modify the numer- 
ical matrix ((26(1 for obtaining a unitary compatible one. 

Taking into account that more often than not unitar- 
ity constraints are not satisfied by experimental data as 
the numerical results ((30|) show we have to devise a fit- 
ting model by supplementing the usual x 2- test with a 
separate piece that should implement the fulfillment of 
unitarity constraints, if we want to get reliable values for 
the interesting physical parameters. 
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III. RECONSTRUCTION OF UNITARY 
MATRICES FROM DATA 

We have seen in the previous section that by using 
numerical values from neutrino data, Eq.Q, one finds 
values for cos S outside the physical range [—1,1], and 
generically all numerical values are different. The ex- 
pressions for cos S are provided by the last four relations 
(|12[) . and these formulas have to give the same number 
when comparing theory with experiment, by supposing 
the data come from a unitary matrix. Their explicit form 
depends on the four independent moduli used to parame- 
terize the data, see e.g. (THI?)) - <j2"T|) . and since there are 58 
such independent groups one gets 165 different expres- 
sions for cos <5. 

Hence if data are compatible to the existence of a uni- 
tary matrix the mixing angles sj^ coming from all 58 

groups, and all 165 phases 5^ m \ expressed as functions 
of V a i, have to be (approximately) equal, and these are 
the most general necessary conditions for unitarity; they 
can be written as 



o<4 m) <i, 



,("0 



» 



cos 5^ = cosS^\ i,j = 1, 



m, n = 1, 
..,165 



,58, (31) 



The above relations are also satisfied by double stochastic 
matrices, as the numerical relations from Eas. (|27[l - il28|) 
show, and the condition that separates unitary matrices 
from double stochastic ones is given by relation i|15[l . i.e. 
-1 < cos<5« < 1. 

As a warning, the content of relations (|31[) can be sum- 
marized as follows: unitarity implies strong correlations 
between all numerical values V a j- 

Now we present the necessary input that have to be 
taken into account when doing a fit the experimental data 
are considered moduli of a unitary matrix. It has to in- 
clude as a separate piece the conditions for a complete 
implementation of unitarity. For that we define a test 
function that has to take into account the double stochas- 
ticity property expressed by the conditions (| 1 31) and the 
fact that in general the numerical values of data are such 
that cos 5 depends on the choice of the four independent 
moduli and could take values outside the physical range. 
The proposal for the first piece is 



X \ = ]T(cos^>-cos^>) 2 + £ £ 

i<j at=e,/x,T \i— 1,2,3 

+ E ( E V 2j-A > -1<cos<5«<1 (32) 

Q=e,/x,r \»=1,2,3 / 

This proposal was made in [3l| for the case of CKM 
quark matrix, and is our contribution to the existing 
methods of reconstructing unitary matrices from exper- 
imental data, when data are compatible with their exis- 
tence; it expresses the full content of unitarity. We stress 



again that both conditions have to be fulfilled: \i nas 
to take small values and cos S values have to be phys- 
ically acceptable, since the relation x\ = holds true 
for both double stochastic and unitary matrices. When 
the experimental data are those from the matrix Q the 
second piece has the form 



(33) 



where V a i — \U a i\ are the minimizing parameters, and 
V a i is the numerical matrix that describes the experimen- 
tal data as in (Q, and a is the errors matrix associated 
to V a i) m our case V e i = 0.835, cr e i = 0.045, and so on. 
Thus the simplest test function has the form 




2 2,2 

X =Xi+X2 



(34) 



Of course unitarity constraints (|32J) could be imple- 
mented by penalty functions, and the concrete form of 
(|34|l could depend on the quality of the fit, which finally 
will choose between the form (|32|) . or that with penalty 
functions. The formula l|34l) provides a least squares 
method which is similar to that used by the BaBar Col- 
laboration j2|| , the difference being the term \\ that im- 
plements the theoretical model (@J. 

The above form can be easily modified if we have other 
experimental information for some functions depending 
on Sij and 5, or on V a i- The primary data in neutrino 
physics are the transition probabilities between two neu- 
trino flavors v a — > i/g, a, (3 = e, //, r at a distance L, 
whose generic form is 



P{y a - v 0] E,L)=J2 U Pl U* P] U al U a3 expH Aj 

= E cos A l + E «n A| (35) 



i>3 



1,3 



where 

W a } = U aj U*^U* m U Pl , A?- = K 2 - m))L/2E v (36) 

and a — e, /j>, v, i,j = 1,2,3. From [l^, see also [33$ . 
we know that imaginary part has the form 

^(UajU^U^Upi) = J ^e a ^e jU (37) 

with J the Jarlskog invariant. On sees that the parame- 
ter involving the CP violation appears in transition prob- 
abilities, hence it can be directly determined from the 
experimental data. By using the above relation one gets 



^(U^U^Upi) sin A 



±2 J(sin A 2 2 + sin A^ + sin A 2 3 ) 



(38) 
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The term coming from the real part ^(W 7 ^) depends where 
explicitly on the flavors a and (3. When a = fi and (3 — e 
one gets 

cobA% = 

i,3 

Rl2 cos A 2 2 + i?i 3 COS Ag X + R 23 COS A23 (39) 



Rl2 = U el U e2 $l(U ftl U* 2 ) = - (c? 2 Ci 3 S? 2 (c| 3 - S? 3 S2 3 ) + Ci2Ci 3 C 23 Sl 2 Sl 3 S23(c? 2 - s\ 2 COS (5) 

#13 = C^el C^e3 C^m3 ^(^Ml) = _C 12 C 13 S 13S23(C12S1 3 S 2 3 + C 23 Si2 COS(5) (40) 

-R23 = U e2 U e3 U^RQJ^) = Ci 3 si 2 si 3 s 23 (ci2C 23 cos(5 - si 2 si 3 s 23 ) 



The above formulae have been obtained by using the 
standard form (gj of the PMNS matrix. We use now 
the phenomenological model i|19fl to write i?y in terms 
of invariant quantities, i.e. moduli. One gets 



R12 = -(-b 2 -a 2 c 2 + b 2 c 2 +e 2 -a 2 e 2 ) 
R 13 = l(b 2 -a 2 c 2 -b 2 c 2 -e 2 + a 2 e 2 ) 
i? 23 = \{-b 2 + a 2 c 2 + b 2 c 2 -e 2 + a 2 e 2 + 2b 2 e 2 ) 



(41) 



Similar formulas have to be obtained for all amplitudes 
entering the experimental data. For exarnple the three- 
flavor survival probability P ee in vacuum Q is given by 

P ee = 1 - i sin 2 26»i3 - cos 4 6 1S sin 2 20 12 sin 2 A 2 21 (42) 
= 1 - 2(a 2 + b 2 ) + 2(a 2 + 6 2 ) 2 - 4 aV sin 2 A 2 i 



Hence the transition probabilities l|4l)|) - (|42[l . and those 
similar to, depend on eight parameters: a, 6, c, e, J, 
A12, Ai 3 , A2 3 which can be used to fit the data. In the 
usual assumption that U is unitary only six are indepen- 
dent, for example a, b, c, J, A12, Ai 3 . That means that 
we may eliminate modulus e. For that we use the rela- 
tion J — 2 A where, J is the Jarlskog invariant, and A is 
the area of (any) unitarity triangle. One gets 

p2 _ b 2 (l - a 2 - c 2 ) + a 2 c 2 (l - a 2 - b 2 ) 

2\2 ' 



(l-« 2 ) 

2ey , a 2 fr 2 c 2 (l - a 2 - c 2 )(l -a 2 - b 2 ) - (1 - a 2 ) 2 J 2 
(1 - a 2 ) 2 



(43) 



where e = ±1. If the matrix is unitary then the following 
condition has to be satisfied < e 2 < 1. The second 
parameter is eliminated from the relation 



If we use the the phenomenological models l|20|) - (|21|l . i.e. 

if we take as parameters four independent moduli, we find 
form different terms for all Rij appearing in formulae 



(gUJ), see equations (gTJ- In- 
form 



Thus x 2 could have the 



X2 



= EEE 



(44) 



A 2 2 + A 2 3 + A 2 j = 



where a runs over all transitions probabilities involved 
in the fit, i runs over the sets of four independent mod- 
uli, such as those given by relations (|19|) - (|21|l . which give 
equivalent but form different expressions for the theoret- 
ical probabilities P(^ Q ,_ >l//3 ), and j runs over the number 
of experimental points for each transition probability. By 
doing a fit with the form (|34|l where x 2 is given by l|44|) 
one gets values for moduli V a i, A^ and J. 

The fit could provide a few (approximate) unitary ma- 
trices compatible with data. In such a case the dou- 
ble stochastic matrices obtained from fit could be used 
to do statistics on unitary matrices. At this point the 
embedding property of unitary matrices into the convex 
set of double stochastic matrices is essential. Indeed if 
Ui, U 2l ... U n are unitary matrices, then 



M 2 = ^xHf/,1 2 , 5> i = l, (45) 

i=l i=l 

< X4 < I, i = 1, . . . , n 



is double stochastic. The above relation shows that the 
statistics on unitary matrices has to be done through the 
statistics on double stochastic matrices, i.e. the relevant 
quantities for doing statistics on unitary matrices are the 
squared moduli, |J7 Q:; | 2 . This property allows us to define 
correctly the mean value, < M >, and the error matrix, 
a M-, for a set of doubly stochastic matrices, coming from 
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a set of approximate unitary matrices, as follows 
<M> 



/k—n 



(46) 



cm 



\ 



S i^'i 4 / n - < M >4 < 47 ) 



\/c=i 



If the moduli of the mean value matrix, < M >, obtained 
in this way are compatible with a unitary matrix, or not 
too far from moduli coming from a unitary one, one can 
reconstruct from it an (approximate) unitary matrix, and 
in the following we use this technique to reconstruct uni- 
tary matrices compatible to the data, and for obtaining 
information on relevant physical parameters. 



requirement see, e.g., Ref. [3^, the aim in view being a 
good estimation of mixing and CP violation in neutrino 
physics. For doing that we use the simple ansatz 
(|34[l in conjunction with data JQ| to obtain qualitative 
results about J, 8 and/or V e3 . The last parameter was 
considered free within the interval V e s G (0.02 — 0.4) and 
its error was taken cr e3 = V e 3 . Around the central values 
from we obtained many approximate double stochas- 
tic matrices, i.e. values for all moduli V a i, and by doing 
statistics with formulae 147fl the following mean value 
matrix and associated error matrix were got 



0.835392 0.537838 0.114933 \ 
Mi = | 0.381844 0.585045 0.715381 
0.395402 0.607056 0.689253 / 



(48) 



IV. GLOBAL FITS 

By global fit we understand a fit that takes into ac- 
count all available experimental data, and (all) unitarity 
constraints expressed by invariant functions; on this last 



0.004241 0.00467 0.009417 
cr Ml = | 0.001008 0.00089 0.000263 
0.003287 0.005701 0.009052 

and the unitary matrix corresponding to Mi is 



(49) 



0.835176 0.537838 0.114933 \ 

Ui= I -0.0697446- 0.375651 i -0.0445706 + 0.583326 i 0.715381 (50) 
-0.0668811 + 0.38991 i -0.0434268 - 0.605467 i 0.68922 / 



If we define two matrices, let say, by the relation 

M± = ^Ml±ial h (51) 

we get the following prediction for CP violating phase 
and unitarity triangle area 

8 = (89.98t[J;26) , A = (127.44 ± 0.01) x 10" 4 

corresponding to a 3 a level, and so on. 

The results (148(1 and (|49jl strongly depend on the mod- 
uli values around one looks for a minimum. Changing 
these values within the error bars in one gets other 
matrices and other values for 8. For example the matrix 

/ 0.851242 0.523442 0.0373568 \ 
M 2 = 0.34007 0.578624 0.741314 (52) 
\ 0.399674 0.625462 0.670118 / 

leads to a value 8 = (120.01 ± 0.16)°, J » 2. x 10~ 4 , 
and so on. This happens since the errors entering 
are too large, such that for many choices of fixed three 



moduli, 8 is practically unconstrained, i.e. there is a 
continuum of unitary matrices with cos 8 varying in many 
cases between (—1, 1). 

The result (|48|l is "confirmed" by a "standard" fit when 
one uses formulas (|12|l , without any approximations com- 
ing from small parameters, etc., and one takes as inde- 
pendent parameters 0ij and <5. The term containing V e 3 
was not introduced in the fit. One gets easily 

0i2 = 0.572332, 6>i3 = 0.088198 , 9 23 = 0.796547, 
8 = 90.51°, A = 9.9x 10~ 3 (53) 

where the values for 9ij are given in radians, and x 2 — 
8.7 x 10~ 2 , where from one gets the following double 
stochastic matrix 

/ 0.837373 0.539489 0.0880842 \ 
M 3 = 0.381844 0.589052 0.712167 (54) 
\ 0.391116 0.601639 0.696462 / 

The corresponding unitary matrix is given by 
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0.837373 0.539489 0.0880842 

U 2 = ( -0.0495463- 0.378656 i -0.0393741 + 0.587735 i 0.712167 | (55) 
-0.0552421 + 0.387195 i -0.0279691 - 0.600988 i 0.696462 



The remarkable fact is that the fit provides a numerical 
value for V e 3 which could not be too far from the true one, 
although the global fit Q failed to find it. To estimate 
the errors to the above parameters we used the moduli 
matrices \JV 2 — 2 a 2 and \JV 2 + 4a 2 , where V and a 
are those given by QJ, and doing again the fit we got 
X2a = 1-32, and, respectively, x\o = 1-22, where from 
the following prediction was obtained 



6= (90.51™)° 



A = (9.2t°;£) x 10~ 3 



(56) 



To see the difference between statistics done by using 
(|47|l and (|51(l and the naive one, we used the central ma- 
trix l[T]l. and the corresponding error matrix a obtaining 
V c + a. to find 



0.57996, 0i3 = 0.00054, 
0.797221, S = 79.54° (57) 

results which are equivalent with the moduli matrix 



712 
923 



0.836484 0.547991 0.00054 \ 
W= ( 0.382939 0.584411 0.715417 
0.391984 0.598473 0.698697 / 



(58) 



This matrix is not compatible with a unitary matrix if 
we use our criteria. For example one gets cos^ 3 -* = 
-0.1046 i, and A = 9.0 x 10~ 5 i. In fact there are 20 cos S 
from a total of 165 which take imaginary values, so, in 
principle, the problem cannot be neglected. On the other 
hand one finds \ 2 — 6-4 for 4 d.o.f., which is not too bad. 
Taking into account that the double stochasticity prop- 
erty (unitarity) is only weakly violated, we find from l|57(l 
a conservative estimate S > 70°. Doing a similar fit for 
V-a one gets \ 2 = 8.1 for 4 d.o.f., and S = 91.81°. 

Of course in true fits one does not work with quanti- 
ties such as V ± a, and they were given as pedagogical 
examples. They show that it is possible when using as 
independent parameters 6ij and d to obtain physical val- 
ues for all of them, although the area of unitarity triangle 
could take imaginary values. The above computations 
show that Sjj and S are convention dependent, even if 
their numerical values are physical, and in order to min- 
imize this bias we expressed these parameters in terms 
of the moduli V a i which have to be the same irrespec- 
tive of the used convention for the unitary matrix 
and we required the fulfillment of all constraints (|31|l . 
The second lesson is that one could get values for 8 even 
if moduli where 5 enters are not yet measured, and this 
happens since unitarity over constrains the data implying 
many relationships between the moduli. Of course a neg- 
ative attitude to such results is compliant so long as the 
subtleties implied by unitarity are not fully understood. 



In the following we want to stress the non uniqueness of 
unitary matrices compatible with the error bars in l[T]l. 
For example, one can find also a unitary matrix which 
lead to 5 = (75.54 ± 0.11)°; its moduli are 



0.813246 0.575234 0.087959 \ 
V I 0.453777 0.614523 0.645328 
0.364304 0.539877 0.758825 / 



(59) 



With the matrices Mi, i = 2,3,4, by using the prop- 
erty of double stochastic matrices we find a continuum 
of (approximate) unitary matrices whose squared moduli 
are given by 

M 2 = x Ml + y M 3 + (1 - x - y) Af 4 2 , < x, y < 1 (60) 

Doing statistics on M,, i = 2,3, 4, with formulae (|47|) one 
gets 

/ 0.834101 0.546484 0.0750359 \ 
M 5 = 0.394701 0.594258 0.700756 (61) 
\ 0.385327 0.590096 0.694424 / 



0.026107 0.023843 0.002995 
om 6 = ( 0.037521 0.0180207 0.05552 
0.01148 0.041863 0.06711 



(62) 



With matrix M one can obtain any value for 8 £ (75° — 
120°), and, respectively, A e (2 x 10~ 4 - 1.6 x 10~ 3 ). The 
prediction at the symmetric point x = y — 1/3 obtained 
from M5 is 



1 n +14.05 

■ i0 -10.45 



)°, A = (11.1 + 0.3) x 10" 



at one sigma level given by <7j\/ 5 . 

In the following we want to test the dependence of 
moduli matrix on the form of the unitary matrix U . For 
that we use the Kobayashi-Maskawa form (2{| 



(Cl C3S1 
C2S1 -C1C2C3 + S2S3 e 1 ' 
S1S2 -C1S2C3 — C2S3 e" 



S1S3 \ 
-C1C2S3 - C3S2 e lv (63) 
-C1S2S3 + C2C3 e 1 ^ J 



where c 2 ; = cos0i and Sj = sin0i, and ip is the Dirac 
phase. Of course selecting one or another form has no 
theoretical significance because all choices are mathemat- 
ically equivalent; however a clever choice may shed some 
light on important issues. Such issues could be the |J/ e 3| 
magnitude and the CP phase 6, (ip) , and the above form 
(|63H could be more helpful. 
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Doing a fit by using the relations (|12f> and the new 
parameters 0j one gets the moduli matrix (|54|l with the 
same digits. In fact we get 

6 l = 0.578337, 2 = 0.797342, (64) 
3 = 0.161845, <^ = 90.15° 

that confirm the previous results, the only change be- 
ing the value for 63 which is brought, in principle, to a 
measurable value. 

A few conclusions of the above numerical fits could be: 
a) there are large corridors for moduli around the central 
values in Q where data are compatible to the existence 
of a unitary matrix with 6 7^ 0, 7r, and by consequence 
compatible with a CP violation; b) the real parts of the 
complex entries in unitary matrices are smaller by an or- 
der of magnitude than the imaginary parts, see e.g. ma- 
trix Ui, but their neglect completely spoils the unitarity 
property; c) by consequence, reliable results can be ob- 
tained then and only then when all theoretical formulas 
used in fits are exact, not approximate. 

We are aware that the "experimental" data JIJ are not 
numbers obtained directly from experiments. In the fol- 
lowing we suggest a new kind of global fit able to resolve 
the existence or the lack of the CP non conservation in 
the leptonic sector by measuring directly the Jarlskog in- 
variant J, and by testing also the possible new physics 
which could explain the negative value for |£/ e 3| 2 obtained 
in 0|, by using the present day experimental data. For do- 
ing this we suggest the use of moduli, a, 6, c, d, e, /, g, h 
and J, A12, A13, as parameters, although not all of them 
arc truly independent. However the moduli are corre- 
lated through different forms for J and/or cos 5, and re- 
lations as (|3ip. or J (a, 6, c, e) = J(c,d,g,h) have to be 
satisfied to a great accuracy. And such a fit could provide 
values for all moduli. If the new fit gives a positive value 
for J one can say that the fit confirms the CP violation 
in the leptonic sector. If the fit will confirm the results 
from Ref. 0, see its Appendix A, concerning the best fit 
obtained for negative values of |C/ e 3| 2 , which imply that 
a and b entering the formula l|42[l for P ee are such that 
1 — a 2 — b 2 < 0, one can use the unitarity constraints 
to obtain other correlations. For example the relation 
(021 implies also that 1 — a 2 — c 2 = \U T i\ 2 < 0, because 
otherwise e 2 gets complex. If this is true then by using 
as independent parameters a, c, e, and/, and eliminating 
again the parameter e, one gets 

2 _ (a 2 + c 2 )(a 2 -c 2 / 2 )-(a 2 ~c 2 ); 2 
(a 2 + c 2 ) 2 

2eVa 2 c 2 / 2 (l - a 2 - c 2 )(a 2 + c 2 - f 2 ) - (a 2 + c 2 ) 2 J 2 

(a 2 + c 2 ) 2 (65) 

Because |kri| 2 = 1 — a 2 — c 2 < 0, the preceding relation 
tell us that |[/ r2 | 2 = a 2 + c 2 - f 2 < implying a new 
physics; in fact a new type of unitarity. In any case by 
using the proposed fit there is a bigger flexibility to test 
various consequences of the theoretical model Q. 

Of course results such as \U Ti \ 2 < 0, i = 1, 2, or 
I t/e3 1 2 < could be artifacts of approximations used for 
the transition probabilities by neglecting the sub leading 



terms, and before speaking of new physics one must use 
exact theoretical forms for all observables measured in 
neutrino physics. The unitarity is a very powerful tool 
because it implies many constraints between the moduli, 
and even a moduli variation of the order of 2-3% could 
spoils it. See in this respect the numerical matrix l|26|l . 
and/or the matrix 1)58(1 . obtained from physical values 
for Sij and <5, which show easily how one can walk into 
a trap. Such a fit will provide values for all moduli even 
if the nowadays data come only from three sources, and 
it deserves to be done since it will give estimates for all 
transition probabilities necessary for the design of future 
neutrino factories. However, if such a result is confirmed 
this will imply a new physics which could be interpreted 
as the first sign showing an evidence for a new generation 
of leptons, or the change of the theoretical frame (0J, by 
taking the matrix U as being unitary, for example, in a 
non Euclidean metric. 



V. CONCLUSIONS 

In this work we have discussed the question of the rel- 
evance of the nowadays experimental neutrino data on 
the determination of CP non-conservation in the leptonic 
sector, as well as on solving the negativity problem of 
I C/ e 3 1 2 < 0, raised by the best fit in Ref. y|. For doing 
that we used an other consequence of unitarity property 
and defined a phcnomcnological model by taking as free 
parameters the moduli of the unitary matrix, which are 
invariant quantities, and which, naturally, leads to the 
separation criterion between the double stochastic ma- 
trices and those arising from unitary ones. This crite- 
rion provided the necessary and sufficient conditions for 
the consistency of experimental data with the theoreti- 
cal model encoded by the PMNS matrix, the strongest 
condition being — 1 < cos<5 < 1. 

We constructed an exact double stochastic matrix, 
which from an experimental point of view is perfectly ac- 
ceptable, but which does not come from a unitary matrix, 
showing the necessity of using the separation criterion in 
phenomenological analyzes in order to obtain physically 
meaningful quantities. 

Taking into account the above result, we provided a re- 
liable algorithm for the reconstruction of a unitary ma- 
trix from experimental data, when they are consistent 
with the theoretical model. The fitting algorithm has 
the form a least squares % 2 -test and have two separate 
pieces: one which implements the unitarity constraints, 
and the other which properly takes into account the er- 
ror affected data. We have also shown how this algorithm 
can be modified to take into account the oscillation prob- 
abilities that are the primary data measured in neutrino 
physics. 

We have used the x 2 -test to obtain information on the 
four parameters entering the standard parametrization 
of the PMNS matrix. The test works very well, and by 
applying it we obtained a continuum of (approximate) 
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unitary matrices compatible with the nowadays available 
experimental data summarized in the numerical matrix 
Q. Their explicit construction shows that the moduli 
data m are compatible with the existence of CP viola- 
tion in the leptonic sector, even if one cannot yet find a 
precise value for phase S. 

The numerical results show that unitarity is a very 
strong property, and the tuning of moduli implied by 
it is given at a higher level than expected, the statistical 
errors generated by the approximate character of data are 
at least an order of magnitude less than the experimental 
errors. 

It was also shown that the fitting method, using as 
independent parameters moduli and Jarlskog invariant, 
can directly check the existence of CP violation in the 
present day neutrino data. On the other hand, the same 
method can resolve the positivity problem for |J7 e 3| 2 . If 
the results by Maltoni et al. [j| are confirmed, they will 
be the first signal for a new physics in leptonic sector. 

Last but not least we used the natural embedding (J7J 
of unitary matrices into the double stochastic set for de- 



vising a method for doing statistics on unitary matrices, 
problem that was an open one until now. The method 
allowed us to find a procedure to obtain mean values, 
and, correspondingly, error matrices starting with a set 
of (approximate) unitary ones. 

As a final conclusion, a satisfactory solution of CP non- 
conservation can be obtained by diminishing the errors 
in all running experiments, as well as in all future ones, 
the present theoretical approach being able to test the 
validity of the nowadays unitary model, or the necessity 
to modify it. 
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